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We propose to measure the CP phase ócp, the magnitude of the neutrino mixing 
matrix element |Ue3| and the sign of the atmopheric scale mass-squared difference 
Am3, with a superbeam by the joint analysis of two different long baseline neutrino 
oscillation experiments. One is a long baseline experiment (LBL) at 300 km and 
the other is a very long baseline (VLBL) experiment at 2100 km. We take the 
neutrino source to be the approved high intensity proton synchrotron, HIPA. The 
neutrino beam for the LBL is the 2-degree off-axis superbeam and for the VLBL, 
a narrow band superbeam. Taking into account all possible errors, we evaluate the 
event rates required and the sensitivities that can be attained for the determination 
of dcp and the sign of Am3,. We arrive at a representative scenario for a reasonably 
precise probe of this part of the neutrino physics. 


I. INTRODUCTION 


The Super-Kamiokande experiments [1] in the past several years, joined by SNO [2] more 
recently, have given strong indications of neutrino oscillation that are corroborated and con- 
strained by a variety of other experiments. These experiments started a new era in the study of 
neutrino physics and offered the best indication to date of physics beyond the standard model. 
To further probe neutrino physics, there are a number of ongoing and planned neutrino oscilla- 
tion experiments. These experiments promise to give a full description of the phenomenology of 
neutrino mixing. The most attractive experiments among the new generation of neutrino oscilla- 
tion experiments are the long baseline (LBL) experiments. They are performed in the controlled 
environment of traditional experimental high energy physics and expected to allow precision mea- 
surements of the oscillation parameters, including the leptonic CP phase. Notably, the recently 
approved superbeam facility, the High Intensity Proton Accelerator (HIPA) [3], which can provide 
intensive high energy neutrino beams from its 50 GeV proton synchrotron, offers the possibility 
of even more desirable LBL experiments. So far, the possibility of two LBL experiments using 
the HIPA superbeam have been discussed. One is HIPA to Kamiokande at a baseline length of 
about 300 km [4] known as J2K, and the other is HIPA to a detector located 2100 km away near 
Beijing [5,6] called H2B. It is well-known that there are parameter ambiguities that are generally 
associated with oscillation measurements at a single baseline [7,8]. Measurements at more than 
one baseline can be beneficial [6,9,10]; our previous studies [11,12] showed that the joint analysis 
of the J2K and H2B experiments can offer extra leverages to resolve some of these ambiguities. 
Our results, however, also showed that CP violation effects cannot be determined at 3o level 
even with the joint analysis considered in the study, in which no antineutrino beams were used. 

Since the leptonic CP phase and mass-squared difference sign are pertinent information in 
the physics of neutrino mixings, which seems to be very different from that of the quark sector, 


it is necessary to find out how to pin down these neutrino mixing parameters accurately in new 
experiments. It has been widely recognized that the neutrino factory [13] is an ideal facility 
for the study of neutrino mixings. However, because of the technical and budgetary challenges 
faced with building a workable neutrino factory in the near future, and because of the availability 
of a conventional superbeam from HIPA in about five years, it is obviously advisable that we 
explore the full potential of the HIPA superbeams. In this work, we examine in further detail the 
measurement of the CP phase and the mass-squared difference sign from the joint analysis at 
the two long baselines with specific HIPA superbeams and more suitable detectors. The neutrino 
beams are the 2-degree off-axis superbeams for the LBL at 300 km and a narrow band superbeam 
for VLBL at 2100 km. We evaluate the event rates and investigate their sensitivities to the CP 
phase and the sign of Am3,. Taking into account all possible experimental errors, we find that 
a fairly precise measurement of the CP phase, the sign of the mass-squared difference and the 
mixing angle 013 is possible but requires: (1) the joint analysis at the two baselines; (2) that 
both a v, beam and a 7, beam are needed at 300 km; (3) that a 7, beam is needed at 2100 km 
if Am}, is negative; and (4) a significant increase in the statistics at both 300 km and 2100 km. 
We find that sin?(20,3) can be probed to very small values, depending on the value of the CP 
phase. 

In Sec. II we describe how the simulations are performed. In Sec. III our results are presented. 
A brief discussion and conclusion can be found in Sec. IV. 


II. DESCRIPTION OF SIMULATIONS 
A. Parametrization and inputs 


Our oscillation analyses will be restricted to 3 flavors of active neutrinos. The parameters of 
the system consists of 2 mass-squared differences (MSD), 3 mixing angles and 1 measurable CP 
phase. The unitary mixing matrix in the vacuum is parameterized as usual 


C12C13 C13512 $13 
U= — €23812 — C12513523  C12C23 — 5125135823 C13823 > (2.1) 
$12523 — C12C23813  —C12823 — C23812813 C13C23 


where sj; = sin(0jx), cjy = cos(0jx), and 3j, = sin(O;,)e°°?, 0;x defined for j < k are the 
mixing angles of mass eigenstates v; and vk, and dcp is the CP phase angle. The three mass 


eigenvalues are denoted as mi, m2, and ma. The two independent MSD are Am3, = m? — m? 


and Am}, = m? — mj. 
The inputs of the mixing angles and MSD’s are obtained from solar, atmospheric and reactor 


experiments: 
sin?(2012) = 0.8, Ami, =5 x 10 ?eV?, 
sin?(2023) = 1.0, |Am3,|=3x1073eV?, (2.3 
sin? (2013) < 0.1 
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Note that the sign of Am2, is unknown. The currently favored Large Mixing Angle solar solution 
requires 6m3, > 0. 

In LBL experiments the neutrino beam interacts with electrons contained in Earth matter [14]. 
In the present study we use the Preliminary Reference Earth Model [15], generally known as 
PREM, for Earth density profiles [16] and numerically integrate the Schródinger equation that 
descibes the propagation of the neutrino in matter for the treatment of distance dependent matter 
density. However, we note that there exist more sophisticated approaches to Earth matter effect, 
including both an updated average density profile known as the AK135 [17] and treatments of 
uncertainties of the density profile [8,18,19] which can affected the determination of the CP phase 
angle. 


B. Beams and detectors 


The HIPA 50 GeV proton synchrotron beam calls for a power of 0.77 MW in phase I, to be 
upgraded to 4 MW in phase II. The superbeam provided by HIPA can be a wide band beam 
(WBB), a pulsed narrow band beam (NBB), or an off-axis beam (OAB). The WBB contains 
neutrinos with widely distributed energy. In a NBB the neutrino flux is concentrated in a narrow 
range of energies, with maximum energy Epeak where the intensity is peaked, and the intensity 
decreases rapidly below Epeax. An OAB also peaks at a certain energy, but has a longer high- 
energy tail than that of a NBB. More details of the various beam profiles can be found in Ref. [20]. 
For the 2100 km baseline we use a NBB [21] with peak energy of 4 GeV. We have also investigated 
the beams of peak energies of 5, 6 and 8 GeV, but found that the 4 GeV beam gives the best 
results. For 300 km we use the 2-degree OAB (2°-OAB) which has a peak energy at 0.8 GeV [22]. 

The detector at 2100 km is assumed to be a water Cerenkov calorimeter with resistive plate 
chambers [23,5] located in Beijing, tentatively called the Beijing Astrophysics and Neutrino 
Detector (BAND). The size of the detector will be 100 kt at the beginning and can be upgraded 
to a much larger one depending on the physics requirements. The detector at 300 km is initially 
the Kamiokande detector of the present size of 22.5 kt and upgraded later to 450 kt.* 


C. Experimental errors 


For the experimental errors we use 30 throughout this work. All independent errors, statistical 
and systematic, are added in quadrature. For the statistical error we used 30 Poisson errors as 
described in the appendix of Ref. [24], including a background at the 1% level of the rate of the 
survival channel. For the systematic error we assumed that the background is known at the 2% 
level as given in [25]. To estimate the error due to the uncertainty in the measurement of the 
mixing angle 023, we assumed that sin? (2023) is measured via the survival channel at L = 300 km, 
with the event rate given by N (v, — v,) ~ No(1 — sin? (2023) sin?(A)), where A ~ 7/2 and No is 
the number of events in the absence of oscillations. Then the statistical uncertainty on sin? (2023) 
is 


§(sin2(203)) = VN/No © = (1 — sin? (2655) sin?(A)). (2.5) 


We then find the variation in the rate in question for a 30 deviation in 055, and added this in 
quadrature to the other 30 errors described above to obtain the total 30 error. 


D. Scenarios 


We consider three scenarios in the present investigation. The scenarios are summarized in 
Table 1. For Scenario I the first stage involves a 5-year experiment with a water Cerenkov 
detector of 22.5 kt at 300 km with the 2°-OAB v,, beam [22] from HIPA at 0.77 MW. This stage 
is contained in the plan of J2K [4]. The second stage of this scenario has HIPA upgraded to 4 
MW, as discussed in Ref. [5], to deliver à NBB v, beam to the water Cerenkov detector of 100 
kt at L=2100 km to run for 5 years. 

Scenario II has an upgraded 4 MW HIPA and calls for both v, and 7, beams. The experiment 
at 2100 km is the same as in Scenario I. For 300 km, however, we assume a much larger water 
Cerenkov detector of 450 kt. It will run for 2 years with 2?—OAB v, beam, and then 6 years 
with the 2°-OAB r,, beam. 


“For a more discussion of the 300 km detector we refer to Ref. [4]. 


Scenario III is similar to Scenario II, but calls for a much larger water Cerenkov detector at 
2100 km, to run either the v, or 7, beam, for example, for 5000 kt-yr. Whether a v or 7, beam 
is delivered to the 2100 km site depends on the sign of Am3,, as will be explained below. 


III. RESULTS OF THE SIMULATION 
A. Strategy 


Let us first describe briefly the strategy of our calculation. There are 6 measurable oscillation 
parameters in the 3-flavor neutrino scheme. There are two mass scales; one is the atmospheric 
scale and the other the solar scale. Existing oscillation experiments have determined that the 
two mass scales are widely separated and therefore sensitive to different L/ E, regions. The LBL 
and VLBL we are considering are affected by the solar scale only at next-to-leading order, but 
can be strongly affected by the unknown parameters dcp and 013. Hence we will take Am}, 
and 012 to be the values determined in solar experiments, as given in Eq. (2.2). This leaves 4 
parameters to be determined, but |Am3,| and 023 are known already to a fair degree of accuracy 
from atmospheric neutrino experiments. Therefore obtaining 013, dcp, and the sign of Am}, is 
the main goal of our calculation. 

The current and soon to be online LBL experiments will determine |Am2,| and sin?(2023) to 
a better accuracy; 6m3, and sin?(2013) will be measured more precisely by KamLAND [26,27]. 
The first task of the scenarios we propose is to determine |Am3,| and sin? (2023) more accurately, 
using the v,, survival probability. This will allow us to have as small an uncertainty as possible in 
the determination of the parameters 013 and dcp. Let us define N,(L) as the number of charge- 
current events involving the a charged lepton at the baseline L. We assume that N,,(300), which 
depends on the v, survival probability, is used to determine |Am3,| and 023 with as small an 
error m possible. Then in the various scenarios N,,(300) cannot be used to determine 013 and 
CP. 


B. Scenario I (v, beam only) 


In this scenario we assume that only the v, beam is employed to run at L=300 km and 
2100 km. Since N,,(300) has already been used to determine |Am3,| and sin?(2623), we are left 
with three types of independent measurements for the determination of 013, dcp and the sign of 
Am3,: Ne(300), Ne(2100) and N,,(2100). The measurements of these three types of events form 
a surface in a three-dimensional space when 013 and óc p are varied in their allowed ranges. The 
angle 013 is constrained by the CHOOZ reactor experiment [28] and the CP phase is completely 
unconstrained. Therefore, we take their ranges to be: sin?(2013) = (0,0.1) and dcp = (0,27). 
Such three-dimensional surfaces, which are tube-like, are displayed in Fig. 1. The upper and 
lower surfaces are for negative and positive Am3,, respectively. The closed curves around the 
axes of the tubes are traced out by varying dcp from 0 to 27, while the lines running parallel to 
the axes of the tubes are determined by varying sin?(2013) from 0.01 to 0.1. For fixed sin? (2013) 
values, we then obtain the ellipses in Fig. 2. 

When N.(300), the number of the ve appearance events at 300 km, is measured, it determines 
a closed curve which is obtained from the three-dimensional surface by a cut at a given value 
on the N.(300) axis. The value of N.(300) does not determine 013 directly since ô is unknown; 
for each of the closed curves we obtain a definite relation between 6 and sin?(2013) when the 


*When new runs at L=300 km with better statistics are made, the improved N, (300) will be used to 
update the values of |Am3,| and 023. 


sign of Am}, is given. We show in Fig. 3 two sets of such relations for each Am3, sign. As 
shown, we choose two extreme values of N.(300), each of which leads to a range of values for 
sin?(2013), depending on the sign of Am?,. For positive Am3,, the larger N. (300) curve limits 
sin?(2013) to the range (0.06, 0.1), while the smaller one corresponds to sin*(20,3) lying in the 
range (0.006, 0.01). Similarly, the ranges of the values of sin?(2013) for negative Am3, can be 
read off from Fig. 3. 

We plot in Fig. 4 the two-dimensional curves with fixed N:(300). Note that the scale of the 
horizontal axis is logarithmic. Ifthe scale was linear, the curves would be ellipses. An open square 
indicates the point on a curve with dcp=0°, solid square 90°, open circle 180°, and solid circle 
270°. We also show a representative 30 error bar for each curve. The assignment of statistical 
and systematic errors has been discussed in the preceding section. The total error is dominated 
mostly by the statistical error. One sees that although the sign of Am3, can be determined at 
the 30 level, there is no sensitivity to the value of the CP phase. In particular, the error in the 
N,,(2100) channel is very large in comparison with the range of variation in the number of events 
when dcp varies; we will encounter similar situation in the next scenario. 


C. Scenario II (v and 7 beams) 


By including the 7, beam aimed at the detector at 300 km, we have two more types of events, 
i.e., Nz(300) and N¿(300). So, in addition to the three-dimensional surface in the N.(300)- 
N.(2100)-N,,(2100) space shown earlier in Fig. 1, we also have surfaces in the spaces N.(300)- 
N=(300)-N7(300) and N«(300)-N2(300)-.N.(2100), as shown in Figs. 5 and 8, respectively. With 
several fixed values of sin?(2013) we obtain the curves shown in Figs. 6 and 9. We plot the 
two-dimensional projections of fixed N.(300)=1000 and 10000 in Figs. 7 and 10. 

We found in Scenario I that the 2100 km data with just a v,, source can determine the sign of 
Am}, at 30, but cannot measure the CP phase (see Fig. 4). In contrast, as shown in Fig. 7, the 
300 km data using both a v and 7 source can determine the CP phase in the ranges (7/2, 37/2) 
or (-7/2,7/2), but cannot distinguish between the two ranges since the measurement is only 
sensitive to sind. Furthermore, unless óc p is close to 7/2 or 37/2 the sign of Am, cannot be 
determined once all of the experimental errors, including the error in the determination of 053, 
are taken into account, leaving a four-fold ambiguity. The problem lies in the fact that Nz (300) is 
used; as already noted in Scenario I, survival data provide poor resolution to the CP phase, and 
the matter effect is small at the relatively short distance of 300 km. Six more three-dimensional 
plots which will contain either N;(300), or N,,(2100), or both, can be made, but they are not 
very useful in the present analysis because they involve the survival data. 

In order to obtain good resolution in the sign of the MSD and to distinguish the two ranges 
of the CP phase as discussed in the preceding section, we have to use data of the electron flavor 
only. Hence we need two experiments with different L/E, ratios and one of them should be a 
VLBL for a good sensitivity to the matter effect. This brings us to the combined analysis of 
Ne(300), N=(300) and N.(2100), as shown in Fig. 10. The sign of Am3, can be easily determined 
if sin2(2013) is not too small and the CP phase can be measured with again the ambiguity 
between the two ranges (7/2,37/2) and (—7/2,7/2), as in the case of Fig. 7. The problem lies 
in the fact that the resolution in N(2100) is poor due to the low number of events, while the 
resolution of the 300 km De is excellent. So we have to increase the statistics at 2100 km. This 
takes us to Scenario III below. 


D. Scenario III (v, and 7, beams with increased statistics) 


The situation of Scenario II can be improved if the statistics in N.(2100) are significantly 
increased. This can be achieved by using a larger detector and/or running for a longer period of 
time for the N.(2100) measurement. For Scenario III we set the detector size times the running 
time at L=2100 km to be 10 times larger than that of Scenario II, assuming that the number of 


events can be straightforwardly scaled up with the detector size. The running at 300 km is the 
same as in Scenario II. The resultant two-dimensional plot is shown Fig. 11. In this scenario, the 
sign of Am3, can be clearly determined at the 30 level, even for N.(300)=1000 which corresponds 
to a very small sin?(20,3) lying in the range (0.006, 0.01), as indicated by the dotted curves in 
Fig. 11. 

If Am, is positive, a reasonably accurate determination of dcp can be made with no 
sgn(Am3,) or 013 ambiguity, and the 7, beam is not needed at 2100 km even for very small 
sin?(2013) in the range (0.006, 0.01). This is consistent with the results of Ref. [9], where it was 
found that a v, — ve and v, — De measurement at short distance and a v,, — Ve measure- 
ment at a long distance could resolve parameter ambiguities for sin?(20,3) > 0.005. To see the 
sensitivity more clearly for positive Am3,, we replot the results in Figs. 12 and 13 respectively 
for N:(300)=10000 and 1000. We see that for N:(300) = 10000, which corresponds to larger 
sin?(2013) (0.06 — 0.1) as shown in Fig. 3, the CP phase can be determined better than 10° at 
30 for dcp small or around 180°. The sensitivity deteriorates slowly when écp moves away from 
0° or 180°, and the uncertainty becomes of the order of 25° when dcp is close to 90° or 270°. 

Even for N.(300)=1000, which corresponds to very small sin? (2013) in the range of (0.006, 0.01), 
the measurement of the CP phase is still reasonably good. It is interesting to note that the 
sensitivity of the CP measurement near écp=0° and 180° for N.(300) = 1000 is comparable to 
that of the much higher number of events of N.(300)=10000. Hence, in this scenario, either case 
can establish whether or not CP in the lepton sector is violated if dcp deviates by than 10° from 
the CP conserving points of dcp=0° or 180°. 

If Am3, is negative, the P, + De oscillation is the favorable channel to investigate. Hence once 
it is clear that Am}; is negative (see the next section for a detailed discussion), the 7,, beam should 
be delivered to 2100 km to run for 5000 kt-yr. The results, which are the counterparts to Fig. 11, 
are shown in Fig. 14. With positron events at 2100 km the CP phase can be well measured. 
To see the sensitivity more clearly, we replot the results in Figs. 15 and 16 for N.(300)=10000 
and 1000 respectively. The accuracy of the dcp measurement for Am3, < 0 using the 7, beam 
is about the same as that of the v, beam for Aime, > 0, although the distinction between dcp 
in the range (-7/2,7/2) and dcp in the range (7/2,37/2) is not as good for a D, beam with 
Ama; < 0. 

We have also done the analysis assuming a v or v NBB of peak energy 5, 6, or 8 GeV is 
delivered to the detector at 2100 km. We found that for these cases the ellipses in Figs. 11 to 
15 are much flatter than for the 4 GeV NBB, so that they do not do as well in resolving the 
degeneracy in dcp. 


IV. CONCLUSION AND DISCUSSION 


We conclude that with a superbeam, such as that delivered by HIPA, the joint analysis at two 
baselines, of which one is an LBL at 300 km and the other a VLBL at 2100 km, can determine 
the Am3, sign and give a reasonably precise measurement of the CP phase and 013. To achieve 
this, both y, and 7, beams are needed for the LBL experiment. The survival events Y, — vy 
and 2, — D, are generally insensitive to the matter and CP effects. 

The initial HIPA v, beam with power 0.77 MW will run with exposure 22.5x5 kt-yr at a 
detector at L—300 km to obtain both survival events N,,(300) and appearance events N.(300). 
The former is used to improve the determination of the mixing angle 023 and mass-squared 
difference |Am3,|, so as to reduce the uncertainty of these crucial input parameters. The latter 
can show the existence of an appearance signal for sin? (2013) > 0.006 and find a crude relation 
between dcp and 013 as shown in Fig. 3. 

A detailed determination of the oscillation parameters will require an upgrade of the HIPA 
beam power to 4 MW. Using our studies in this paper as a guide, we suggest as one possibility 
the following experimental steps using the upgraded HIPA beam: 


Stage 1: Deliver a 4 MW v 2?-OAB to a 450 kt detector at a distance of L=300 km for 2 
years. The survival events N,,(300) are used to determine more precisely the parameters 623 and 


|Am3,|. The appearance events N.(300) are used to refine the relation between dcp and 013, as 
shown in Fig. 3. 

Stage 2: A 4 MW y NBB with peak energy around 4 GeV is delivered to a detector at L=2100 
km, to run for 100x5 kt-yr. The survival and appearance events N,,(2100) and N.(2100) are 
used to determine the sign of Am3,, with the most sensitivity coming from N.(2100) (see Fig. 4). 

Stage 3: A 4 MW D, 2°-OAB is delivered to the 300 km baseline detector for 450x6 kt-yr 
and N¿(300) and Nz(300) are obtained. The data can only determine dcp and 0 up to a 2-fold 
degeneracy because of the poor separation between 6 and 7 — ô in the Nz; (300) measurement, as 
demonstrated in Fig. 7. 

Stage 4: A 4 MW v, (7,) NBB with peak energy around 4 GeV is delivered to the 2100 km 
baseline detector for 1000x5 kt-yr if Am3, > 0 (Am3, < 0). Then at 30, the value of dcp can 
be determined to about 10° for values close to 0° or 180°, or to about 25° for values close to 
90° or 270°. The distinction between dcp in the ranges (-7/2,7/2) and (7/2,37/2) is better 
for Am3, > 0 than for Am3, < 0. The Am}, > 0 case is shown in Figs. 12 and 13 and that of 
Am, < 0 in Figs. 15 and 16. 


It is apparent from our calculation that in order to obtain enough statistics to provide a 
reasonably precise measurement of 013 and dcp the total detector size and running time have 
to be sufficiently large. We have not attempted a detailed optimization; rather, we offer our 
calculation as an example for illustration. A search is still required to determine the optimal 
conditions for the measurement. Eventually uncertainties in the Earth matter density along a 
given baseline as well as uncertainties in the solar neutrino oscillation parameters 012 and Am}, 
must also be taken into account. 
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Table 1 Different possible scenarios in joint analyses 


beam  |power|detector size x runing time] beam | power [detector size x iù time 
(2°-OAB)|(MW) (kt- year) (NBB)|(MW) 
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FIG. 1. Three-dimensional surface in the events space Ne(300) — N-(2100) — N,(2100) in Scenario 
II with CP phase dcp varying from 0 to 27 and sin?(2013) from 0.01 to 0.1. The lower (upper) one 


is for Am32 > 0 (< 0). The surface in Scenario I is obatined by scaling Ne(300) axis by a factor 
0.77(MW) x 5(year) X22.5(kt) „~ 1. 
A(MW) x 2(year)x450(kt) ~ 43° 
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FIG. 2. Same as Fig. 1, but for fixed sin?(2013) = 0.02, 0.04, 0.06, 0.08 and 0.1 for the elipses from 
left to right. 
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FIG. 3. CP phase dcp (in degrees) versus sin?(2013) for fixed N-(300). 
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FIG. 4. N,(2100) versus Ne(2100) for fixed Ne(300). The open square, open circle, filled square and 
filled circle denote dcp = 0, 7, 7/2 and 37/2, respectively. 30 error bars at some points are also plotted. 
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FIG. 5. Three-dimensional surface in the events space N.(300) — Nz(300) — Na (300) in Scenario II 
with CP phase dcp varying from 0 to 27 and sin?(2013) from 0.01 to 0.1. The upper (lower) one is for 
Amz > 0 (< 0). 


14 


5000 
10000 
15000 
168500 ^ 
| SOA 
168000 
167500 
Ng ( 300km) 
167000 \ Q N 
166500 
15000 
10000 
5000 
N- (300km) 


FIG. 6. Same as Fig. 5, but for fixed sin? (2813) = 0.02, 0.04, 0.06, 0.08, 0.1 for the elipses from right 
to left. 


15 


N, (300km) 


x 10? 


1640 . 
| Scenario || | 
| Ne(300km) = 10000 | 
L MOTION ld | 
| Cee "7 | 

1630 — Am?,,>0 | i 
| t Am 31 >0 : 

1620 — | 
EET — ae: 

AT aU E | 

1610 — | : 

1600 MENTES uu 

10? » 


\:(300km) 


FIG. 7. Nx(300) versus Ne(300) for fixed N.(300) in Scenario II. The open square, open circle, filled 
square and filled circle denote dcp = 0, 7, 7/2 and 37/2, respectively. 30 error bars at some points are 
also plotted. The dashed lines denote errors caused by 023 uncertainty. Note that only the upper error 
bar in Ng changed since we take 023 = 7/4 and any deviation from this always moves the Nj result in 
the same direction; this is an artifact of choosing maximal mixing as our starting point. 
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FIG. 8. Three-dimensional surface in the events space N.(300) — N-(2100) — N2(300) in Scenario II 
with CP phase ócp varying from 0 to 27 and sin?(2013) from 0.01 to 0.1. The right (left) one is for 
Amz > 0 (< 0). 
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FIG. 9. Same as Fig. 8, but for fixed sin?(2013) = 0.01, 0.02, 0.03, - - -, 0.1 for the elipses from the 


lower-right corner to the upper-left corner. 
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FIG. 10. N2(300) versus N«(2100) for fixed N.(300) in Scenario II. The open square, open circle, filled 
square and filled circle denote dcp = 0, 7, 7/2 and 37/2, respectively. 30 error bars at some points are 
also plotted. 
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FIG. 11. Same as Fig. 10, but for Scenario III. 
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FIG. 12. Same as Fig. 11, but for N:(300) = 10000 and 6m3, > 0. The stars denote 10° step of CP 
phase from 0 to 7/2. 
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FIG. 13. Same as Fig. 11, but for Ne(300) = 1000 and ómà; > 0. The stars denote 10° step of CP 
phase from 0 to 7/2. 
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FIG. 14. Same as Fig. 11, but for Nz(300) versus Ne(2100). 
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FIG. 15. Same as Fig. 14, but for Ne(300) = 10000 and 6m3, < 0. The stars denote 10° step of CP 
phase from 7/2 to r. 


24 


3000 
| Scenario Il 
| Ne(300km)=1000 
2500 — 
| Pa 
> 2000 | x” 
xX - X 
O 5 P 4 
e - E 
4500 - ra 
lw L x 
i E 
1000 — e 9 
L P si 
| S 
500 n D p 
| eo 
0 E | | | | | | | | | | | | | | | | | | | | | | | | | | | d | | | | | 
0 50 100 150 200 250 300 350 
\.(2100km) 


FIG. 16. Same as Fig. 14, but for Ne(300) = 1000 and óm2, < 0. The stars denote 10° step of CP 
phase from 7/2 to r. 
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